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Abstract.  In  this  work,  we  first  prove  the  existence  and  uniqueness  of  a  strong 
solution  to  stochastic  GOY  model  of  turbulence  with  a  small  multiplicative 
noise.  Then  using  the  weak  convergence  approach,  Laplace  principle  for  solu¬ 
tions  of  the  stochastic  GOY  model  is  established  in  certain  Polish  space. 
Thus  a  Wentzell-Freidlin  type  large  deviation  principle  is  established  utilizing 
certain  results  by  Varadhan  and  Bryc. 
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1.  Introduction 

The  large  deviations  theory  is  among  the  most  classical  areas  in  probability  theory 
with  many  deep  developments  and  applications  [10,11,13,17,35,37].  Although  it 
appears  to  be  no  literature  on  the  large  deviation  results  for  stochastic  shell  model 
of  turbulence,  a  few  authors  have  proved  the  Wentzell-Freidlin  type  large  devi¬ 
ations  for  the  two-dimensional  stochastic  Navier-Stokes  equations  with  additive 
noise  (e.g.  [6])  and  also  for  multiplicative  noise  (e.g.  [34]).  For  Donsker-Varadhan 
type  large  deviation  study  related  to  Navier-Stokes  equations  we  refer  the  readers 
to  Quastel  and  Yau  [32] .  For  the  treatment  related  to  stochastic  two-dimensional 
vorticity  equations  see  the  work  of  Amirdjanova  and  Xiong  [1],  Several  authors 
have  established  the  Wentzell-Freidlin  type  large  deviation  estimates  for  a  class  of 
infinite  dimensional  stochastic  differential  equations  (see  for  e.g.,  [4,7,9,20,33]). 
In  these  works  the  proofs  of  large  deviation  principle  (LDP)  (see  Definition  2.2) 
usually  rely  on  first  approximating  the  original  problem  by  time-discretization  so 
that  LDP  can  be  shown  for  the  resulting  simpler  problems  via  contraction  prin¬ 
ciple,  and  then  showing  that  LDP  holds  in  the  limit.  The  discretization  method 
to  establish  LDP  was  introduced  by  Wentzell  and  Freidlin[17]. 
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Dupuis  and  Ellis  [13]  have  combined  weak  convergence  methods  to  the 
stochastic  control  approach  developed  earlier  by  Fleming  [16]  to  the  large  devia¬ 
tions  theory.  Our  work  is  based  on  the  theory  introduced  by  Budhiraja  and  Dupuis 
[4] ,  where  they  used  the  stochastic  control  and  weak  convergence  approach  to  ob¬ 
tain  the  LDP  for  the  family  {C/e(W(-))}e>o,  where  Qe  is  an  appropriate  family  of 
measurable  maps  from  the  Wiener  space  to  some  Polish  space.  Their  work  relied 
on  the  fact  that  the  LDP  is  equivalent  to  Laplace  principle  (see  Definition  2.6) 
if  the  underlying  space  is  Polish,  which  is  in  turn  a  consequence  of  Varadhan’s 
Lemma  (see  Lemma  2.7)  and  Bryc’s  converse  to  Varadhan’s  Lemma  (see  Lemma 
2.8).  We  refer  the  reader  to  Theorems  1.2.1  and  1.2.3  in  Dupuis  and  Ellis  [13]. 

This  work  deals  with  an  infinite  dimensional  shell  model,  a  mathematical 
turbulence  model  that  received  increasing  attention  in  recent  years.  Apparently 
there  are  only  two  previous  rigorous  works  on  infinite  dimensional  shell  model, 
namely  [8]  and  [2],  one  in  the  deterministic  case  and  the  other  in  the  stochas¬ 
tic  case  with  additive  noise,  respectively.  In  both  of  these  works  a  variational 
semigroup  formulation  has  been  introduced.  Our  present  work  deals  with  a  more 
general  stochastic  model,  with  multiplicative  noise:  the  proofs  of  existence  and 
uniqueness  of  strong  solutions  are  considerably  more  difficult  in  this  case.  The 
main  result  of  this  paper  is  to  prove  a  large  deviation  principle  for  the  solution  of 
the  shell  model. 

The  construction  of  the  paper  is  as  follows.  In  the  next  section,  we  give  some 
definitions  and  basic  properties  from  the  large  deviation  theory.  In  later  part  of 
this  section,  we  describe  briefly  the  work  of  Budhiraja  and  Dupuis  [4]  to  set  up 
the  ground  for  our  main  work.  In  Sect.  3,  we  formulate  the  abstract  stochastic 
GOY  model  when  the  noise  coefficient  is  small.  We  then  prove  certain  a  priori 
energy  estimates  with  exponential  weight.  These  estimates  together  with  the  local 
monotonicity  property  of  the  sum  of  the  linear  and  non  linear  operators  play  a 
fundamental  role  to  prove  the  existence  and  uniqueness  of  the  strong  solution.  In 
the  last  section,  we  establish  the  LDP  for  the  stochastic  GOY  model  perturbed 
by  a  small  multiplicative  noise. 

2.  Large  deviation  principle 

In  this  section,  we  will  give  an  abstract  formulation  and  basic  properties  for  a 
class  of  large  deviation  problems.  Let  us  denote  by  X  a  complete  separable  met¬ 
ric  space  and  {PE  :  e  >  0}  a  family  of  probability  measures  on  the  Borel  subsets 
of  X. 

Definition  2.1.  A  function  I  :  X  — >  [0,  oo]  is  called  a  rate  function  if  /  is  lower 
semicontinuous.  A  rate  function  I  is  called  a  good  rate  function  if  for  arbitrary 
M  e  [0,oo),  the  level  set  Km  =  {x  :  I(x)  <  M}  is  compact  in  X 

Definition  2.2.  (Large  Deviation  Principle)  We  say  that  a  family  of  probability 
measures  {PE  :  e  >  0}  satisfies  the  large  deviation  principle  (LDP)  with  a  good 
rate  function  I  satisfying, 


Vol.  16  (2009)  Large  deviations  for  the  stochastic  shell  model  of  turbulence 


495 


(i)  for  each  closed  set  F  C  X 


limsup  e  log  Pe{F)  <  —  inf  I(x), 

£—>0  X^F 

(ii)  for  each  open  set  G  C  X 

liminf  e  log  Pe{G)  >  —  inf  I{x). 

£ — ^0  x£G 

Remark  2.3.  For  any  given  {Pe  :  e  >  0}  there  is  at  most  one  rate  function  gov¬ 
erning  the  large  deviations  of  {Pe  :  e  >  0}. 

Example  2.4.  As  a  simplest  example  we  recall  the  one  dimensional  version  of 
Cramer’s  theorem  [37].  Let  Xn  =  (Yi  +  Y2  H - b  Yn)/n  where  {Yj}  are  indepen¬ 

dent  random  variables  with  a  common  distribution  a.  Assume  that  the  moment 
generating  function 

M(0)  =  E[exp(9Y)]  =  j  e6yd a{y) 

is  finite  for  all  9.  Then  {A !„}  satisfies  the  LDP  with  rate  function  (see  [11]) 

I  (x)  =  sup[0x  —  log  M  (0)] . 
e 

Example  2.5.  We  choose  Schilder’s  theorem  as  second  example  which  has  many 
important  applications  in  large  deviation  theory,  such  as,  in  the  derivation  of  the 
Strassen’s  renowned  Law  of  Iterated  Logarithm,  in  the  Wentzell  and  Freidlin’s 
estimate  on  the  large  deviations  of  randomly  perturbed  dynamical  systems,  to 
name  a  few. 

Let  d  £  Z+  and 


A0  = 

For  <p  £  Aq  define 


\m 


4>  £  C([0,  oo);  Rd)  :  </>(0)  =  0  and  lim 

£->-oo  t 


lull  1^)1 

ll^lUo  =  SUp——. 

t>  0  1  +  t 


=  0 


Then  notice  that  Aq  is  a  separable  real  Banach  space  [13].  Next,  we  define 
H1  =  H1([0, oo); Kd)  to  be  the  space  of  4>  £  A0  with  the  property  that  cj){t)  = 
/o  </>(s)ds,t  >  0,  for  some  ^  £  L2([0,oo);Kd)  and  set  ||</>||Hi  =  ||^||l=([0,oo);K“).  for 
<j)£  H1. 

Now  for  given  T  >  0,  we  define  It  ■  A0  — >  [0,  oo]  by 

/t(a)  =  /  i  Jo  \m\2dt  if  V  e  H1, 
oo  otherwise. 

Let  {W(t)}  be  a  standard  Wiener  process  in  Rd.  Let  the  process 


Xn(t)  = 

vn 


1 
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takes  values  in  a  Polish  space  E.  Then  {Xn}  satisfies  the  LDP  on  E  with  the  rate 
function  It  (see  [13]). 


Definition  2.6.  (Laplace  Principle)  For  h  £  C&([0, 1]), 

1  f1 

lim  —log  /  e~nh(D  ds  =  —  min  h(x).  (2-1) 

n^oon  J o  xG  [0,1] 


Lemma  2.7.  (Varadhan’s  Lemma  [37])  Let  E  be  a  Polish  space  and  {Xe  :  e  >  0}  be 
a  family  of  E -valued  random  elements  satisfying  LDP  with  rate  function  I.  Then 
[Xe  :  e  >  01  satisfies  the  Laplace  principle  on  E  with  the  same  rate  function  I  if 
for  all  h  £  Cb(E), 


lim  e  log  E 

e^O 


inf  {h(x)  +  I(x)}.  (2.2) 

xGE 


Lemma  2.8.  (Bryc’s  Lemma  [13])  The  Laplace  principle  implies  the  LDP  with  the 
same  rate  function.  More  precisely,  if  {X£  :  e  >  0}  satisfies  the  Laplace  principle 
on  the  Polish  space  E  with  the  rate  function  I  and  if  the  limit 


lim  e  log  E 

e^O 


-  inf- {h{x)  +/(x)} 

xGE 


is  valid  for  all  h  £  Cb(E),  then  {Xs  :  e  >  0}  satisfies  the  LDP  on  E  with  rate 
function  I . 


Note  that,  Varadhan’s  Lemma  together  with  Bryc’s  converse  of  Varadhan’s 
Lemma  state  that  for  Polish  space  valued  random  elements,  the  Laplace  principle 
and  the  large  deviation  principle  are  equivalent. 

Let  (12,  T ,  P)  be  a  probability  space  equipped  with  an  increasing  family 
{Pt}o<t<T  of  sub-sigma-fields  of  T  satisfying  the  usual  conditions  of  right  conti¬ 
nuity  and  P-completeness.  Let  H  be  a  real  separable  Hilbert  space  and  Q  be  a 
strictly  positive,  symmetric,  trace  class  operator  on  H. 

Definition  2.9.  A  stochastic  process  {W(t)} 0<t<T  is  said  to  be  an  iL-valued  Tt- 
adapted  Wiener  process  with  covariance  operator  Q  if 

(1)  For  each  non-zero  h  £  H,  \Q1t2h\~1(W{t),  h)  is  a  standard  one-dimensional 
Wiener  process, 

(2)  For  any  h  £  H,  (W(t),  h)  is  a  martingale  adapted  to  Tt. 


If  W  is  a  an  LZ-valued  Wiener  process  with  covariance  operator  Q  with 
Tr  Q  <  oo,  then  W  is  a  Gaussian  process  on  H  and 

E{W{t))  =  0,  Co v(W(t))  =  tQ,  t>  0. 

Let  Hq  =  Q 1/2H.  Then  Ho  is  a  Hilbert  space  equipped  with  the  inner  product 
(•>  Oo, 

(u,  v)0  =  {Q~1/2u,Q~1/2v),  \/u,v  £  H0, 

where  Q _1/2  is  the  pseudo-inverse  of  Q 1f2.  Since  Q  is  a  trace  class  operator,  the 
imbedding  of  Hq  in  H  is  Hilbert-Schmidt. 
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Let  Lq  denote  the  space  of  linear  operators  S  such  that  SQ 1/2  is  a  Hilbert- 
Schmidt  operator  from  H  to  H.  Define  the  norm  on  the  space  L q  by  |5|2Q  = 
Tr  (SQS*). 

Let 

A  =  |iL0-valued  {ft}  —  predict  ableprocesses  v  such  that  J  |u(s)|gds  <  oo  a.s|  . 
Define  the  set  Sn  of  bounded  deterministic  controls  as, 


Sn=  |ueL2([0,T];7Jo)  :  jf  |«(a)| 


gds  <  N 


The  set  Sn  endowed  with  the  weak  topology  on  L2([0,  T\;  H0)  is  a  Polish  space 

[12]- 

Define  An  as  the  set  of  bounded  stochastic  controls  by 
An  =  {v  G  A  :  v(uj)  G  <Sjv,P-a.s.}  . 

Let  E  denote  a  Polish  space,  and  for  e  >  0  let  Q£  :  C([0,T];iL)  — >  E  be  a 
measurable  map.  Define 

xe  =  gE(w(-)). 

We  are  interested  in  the  large  deviation  principle  for  Xe  as  e  — >  0. 

Assumption  2.10.  There  exists  a  measurable  map  Q°  :  C([0,T]  :  H)  — >  E  such 
that  the  following  hold: 

1.  Let  {ve  :  e  >  0}  C  Am  for  some  M  <  oo.  Let  Vs  converge  in  distribution  as 
SW-valued  random  elements  to  v.  Then  ge{W{-)  +  fQ  ve(s)ds)  converges 
in  distribution  to  G°(f0  v(s)ds). 

2.  For  every  M  <  oo,  the  set 


u(s)ds  Sm 


Km  =  |S° 

is  a  compact  subset  of  E. 
For  each  /  G  E,  define 


Hf)  =  inf 

j>GL2([0 ,T]-,H0y.f=G°(I0  v(s)dS)} 


{U 


)lods 


(2.3) 


where  infnnum  over  an  empty  set  is  taken  as  oo. 

We  now  state  an  important  result  by  Budhiraja  and  Dupuis  [4], 

Theorem  2.11.  Let  Xe  =  QS(W(-)).  If  {Ge}  satisfies  the  Assumption  (2.10),  then 
the  family  {XE  :  e  >  0}  satisfies  the  Laplace  principle  in  E  with  rate  function  I 
given  by  (2.3). 

Remark  2.12.  1.  Notice  that  since  the  underlying  space  E  is  Polish,  the  family 

{ Xe  :  e  >  0}  satisfies  the  LDP  in  E  with  the  same  rate  function  I. 


498 


U.  Manna  et  al. 


NoDEA 


2.  Assumption  1.  is  a  statement  on  the  weak  convergence  of  a  certain  family 
of  random  variables  and  is  at  the  core  of  weak  convergence  approach  to  the 
study  of  large  deviations.  Assumption  2.  essentially  says  that  the  level  sets 
of  the  rate  function  are  compact. 


3.  The  stochastic  GOY  model  of  turbulence 

The  GOY  model  (Gledger-Ohkitani-Yamada)  [29]  is  a  particular  case  of  so  called 
“Shell  model”  (see  [18]).  This  model  is  the  Navier-Stokes  equation  written  in  the 
Fourier  space  where  the  interaction  between  different  modes  is  preserved  between 
nearest  modes.  To  be  precise,  the  GOY  model  describes  a  one-dimensional  cas¬ 
cade  of  energies  among  an  infinite  sequence  of  complex  velocities,  {wn(t)},  on  a 
one-dimensional  sequence  of  wave  numbers 

kn  =  kQ2n,  k0  >  0,  n=l,2, ... 

where  the  discrete  index  n  is  referred  to  as  the  “shell  index”.  The  equations  of 
motion  of  the  stochastic  GOY  model  of  turbulence  have  the  form 

+  vk2nun  +  i  (aknu*+1u*+2  +  bkn-w^u*^ 

+ckn-2Ukn_lukn_2)  =  fn  +  a n(t,un)^U^  ,  for  71  =  1,  2,  .  .  .  ,  (3.1) 

along  with  the  boundary  conditions 

u_  i  =  Wo  =  0.  (3-2) 

Here  w*  denotes  the  complex  conjugate  of  un,  v  >  0  is  the  kinematic  viscosity 
and  fn  is  the  Fourier  component  of  the  forcing,  a,  b  and  c  are  real  parameters 
such  that  energy  conservation  condition  a  +  b  +  c  =  0  holds  (see  [19,29]).  For  the 
standard  model  a  =  —1,  b  =  1/2  and  c  =  1/2.  For  each  n,  wn  is  one-dimensional 
Brownian  motion  and  the  noise  coefficient  er„  is  assumed  to  satisfy  the  following 
properties, 

a.l.  For  all  t  £  [0,T],  there  exists  a  positive  constant  A'i  independent  of  t  and  n 
such  that, 

Wn(t,un)\2  <  Kikl\un\2, 

a. 2.  For  all  t  £  [0,T],  there  exists  a  positive  constant  A'2  independent  of  t  and  n 
such  that, 

|  (^;Tin)  ^  ni^  1  ^  n)\“  A  2  kn  \  un  Vn  |  . 

3.1.  Functional  setting 

Let  H  be  a  real  Hilbert  space  such  that 

OO 

u  =  (Ml,  u2,  •  ■  ■)  £  C°°  :  ^  | wn | 2  <  00 

n— 1 


H  := 
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For  every  u,v  £  H,  the  scalar  product  (•,  •)  and  norm  |  •  |  are  defined  on  H  as 

OO  /  OC  \ 


(u,v)H  =  Re  H=(^|Mn|2 

n—1  \n= 1 

Let  us  now  define  the  space 

V  :=  |  u  £  H  :  ^  k2n\un\ 2  <  oo  i  , 


which  is  a  Hilbert  space  equipped  with  the  norm 


=  £* 


1/2 


■2WI2 


\n—  1 


The  linear  operator  A  :  D(A)  — >  H  is  a  positive  definite,  self  adjoint  linear 
operator  defined  by 

Au  =  ((du)i,  (Au)2,  ■  ■ ■  ),  where  (Au)n  =  fc2un,  \/uGD(A).  (3.3) 

The  domain  of  A,  D(A)  C  H,  is  a  Hilbert  space  equipped  with  the  norm 

/  oo  \  1/2 

\Md(a)  =  \Au\  =  \y^k*\un\2  j  ,  \/u£D(A). 


Since  the  operator  A  is  positive  definite,  we  can  define  the  power  A1/2  , 

Alt2u  =  (fcitti,  k2u2,  ■  ■  •),  Vm  =  (ui,U2,  ■  ■  •)• 

Furthermore,  we  define  the  space 

D{A1/2)  =  |u=  (iii,u2,...)  :  X>«W2  <  00 1 

which  is  a  Hilbert  space  equipped  with  the  scalar  product 

(u,v)D(Ai/2)  =  (A1/2u,A1/2v),  Vu,  v  £  D(A1/2), 

and  the  norm 

/  °o  \  1/2 

IMb(AV=)  =  (  Ylkn\Un\2  J 

Note  that  V  =  D^A1/2).  We  consider  V'  =  D(H-1/2)  as  the  dual  space  of  V. 
Then  the  following  inclusion  holds 


V  C  H  =  H'  c  V'. 
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We  will  now  introduce  the  sequence  spaces  analogue  to  Sobolev  functional  spaces. 
For  1  <  p  <  oo  and  s£l 


i/p 


Ws’p  := 
and  for  p  =  oo 


=  («i,  «2,  •  •  ■)  :  \\As/2u\\p  =  )"  <  00  h 


\n—l 


Ws’°°  :=  {w  =  (rti,  w2,  ■  ■  :  ms/2w||oo  =  sup  {ksn\un\)  <  oo}  , 

l  l<n<oo  J 


where  for  u  £  Ws,p  the  norm  is  defined  as 


IMIw-  =  \\As/2u\\p. 

Here  ||  •  ||  denotes  the  usual  norm  in  the  lp  sequence  space.  It  is  clear  from  the 
above  definitions  that  W1,2  =  V  =  D(A1/l2). 


Remark  3.1.  For  the  shell  model  we  can  reasonably  assume  that  the  complex 
velocities  un  are  such  that  \un\  <  1  for  almost  all  n.  Then 

OO  /  OO 

INI  %  =  lM™l4  -  ( 

n—1  \n— 1 

which  leads  to  H  C  /4. 


We  now  state  a  Lemma  which  is  useful  in  this  work.  We  omit  the  proof  since 
it  is  quite  simple. 


Lemma  3.2.  For  any  smooth  function  u  €  H ,  the  following  holds: 

\\u\\%  <  C\u\ 2  IMI2.  (3.4) 


3.2.  Properties  of  the  linear  and  nonlinear  operators 

We  define  the  bilinear  operator  :  V  x  H  — >  H  as 

B(u,v)  =  (Bi(u,v),B2(u,v),...), 


where 


Bn{u,V )  =  ikn  (  -  ^«+W*+2  +  <+2<+i)  +  ■ 


In  other  words,  if  be  a  orthonormal  basis  of  H,  i.e.  all  the  entries  of  en 

are  zero  except  at  the  place  n  it  is  equal  to  1,  then 


B(u,v)  =  i  Y  kn  (  -u. 

n—1  ^ 


n+l^n-l  .,(.un+lvn+2 


^n+2^n+l 


l)  +  oUn-lVn~: 


(3.5) 


The  following  lemma  says  that  B[u,v)  makes  sense  as  an  element  of  H, 
whenever  u  €  V  and  v  £  H  or  u  £  H  and  v  £  V.  It  also  says  that  B  (u,  v)  makes 
sense  as  an  element  of  V' .  Here  we  state  the  following  lemma  which  has  been 
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proved  in  Constantin  et  al.  [8]  for  the  Sabra  shell  model,  but  one  can  also  prove 
the  similar  estimates  for  the  GOY  model  (see  [2]). 


Lemma  3.3.  (i)  There  exist  constants  C\  >  0,62  >  0, 


\B(u,v)\  <  Ci  u|||u  , 

Vu  £  V,  v  £  H, 

(3.6) 

and 

\B(u,v)\  <  C2MIMI, 

Vu  £  H,  v  £  V. 

(3.7) 

(ii) 

B  :  H  x  H  — >  V'  is  a  bounded  bilinear  operator  and  for  a 

constant  C3  >  0 

\\B(u,v)\\v  <  C3\u\\v 

|,  Vu,  v  £  H. 

(3.8) 

(iii) 

B  :  H  x  D{A)  — >■  V  is  a  bounded  bilinear  operator  and  for 

a  constant  C4  >  0 

\\B{u,v)\\v  <  C4\u\\Av\, 

Vu  £  H,v  £  D(A). 

(3.9) 

(iv) 

For  every  u  £  V  and  v  £  H 

(B(u,v),v) 

=  0. 

(3.10) 

We  now  present  one  more  important  property  of  the  nonlinear  operator  B 
in  the  following  lemma  which  will  play  important  role  in  the  later  part  of  this 
section.  The  proof  is  straightforward  and  uses  the  bilinearity  property  of  B. 


Lemma  3.4.  If  w  =  u  —  v,  then 

B(u,  u)  —  B(v,  v)  =  B(v,  w)  +  B(w,  v)  +  B{w,  w). 

With  above  functional  setting  and  following  the  classical  treatment  of  the 
Navier-Stokes  equation,  and  in  order  to  simplify  the  notation  one  can  write  the 
stochastic  GOY  model  of  turbulence  (3.1)  in  a  Hilbert  space  H  in  the  following 
way, 

du  +  [vAu  +  B(u ,  u)]  d t  =  f(t)dt  +  aft,  u)dW ( t )  (3-11) 

«(0)  =  uo,  (3.12) 

where  u  =  {u\,u2, . . .)  £  H ,  the  operators  A  and  B  are  defined  through  (3.3) 
and  (3.5),  respectively,  /  =  (fi,  f2, . .  .),a(t,u)  =  (<Ti(t,  uf),  a2(t,  u2),  ■  •  ■),  and 
W  =  (wi,w2, . . .).  Here  (W(f)t>o)  is  a  H- valued  Wiener  process  with  trace  class 
covariance,  and  the  space  Lq  has  been  defined  in  Sect.  2.  The  noise  coefficient 
a  :  [0,  T]  x  V  — >  Lq(Hq-:  H)  is  such  that  it  satisfies  the  following  hypotheses: 

A.l.  The  function  cr  <E  C([0,T]  x  V;Lq(H0;  H)) 

A. 2.  For  all  t  £  (0 ,T),  there  exists  a  positive  constant  K  such  that  \<j(t,u)\2L  < 

K{  l  +  IMI2). 

A. 3.  For  all  t  £  (0,  T),  there  exists  a  positive  constant  L  such  that  for  all  u,  v  £  V, 
W (t,u)  -cr(t,v)\lQ  <  L\\u  —  u||2. 
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Remark  3.5.  The  above  hypotheses  can  be  verified  from  the  assumptions  (A  1  — 
A. 2)  on  the  noise  coefficients.  Notice  that,  Q  :  H  — >  H  is  a  trace  class  covariance 
(nuclear)  operator  and  hence  compact.  So  Hq  =  QX!2H  is  a  separable  Hilbert 
space  and  the  imbedding  of  H0  in  H  is  Hilbert-Schmidt.  Let  { en }^=1  be  the 
eigenfunctions  of  Q  (may  not  be  complete).  Then  Qen  =  \nen,  where  each  Ara  is 
positive  real  and  A„  <  oo.  Note, 

OO  OO 

|cr(£, u)| Lq  =  'y  ~  | (o’/imj £n) I  =  y  '  Am|(crem, en)\ 

m,n— 1  m,n=l 

=  (crQl/2  ,aQ1/2)  =  Tr(aQa), 

where  {hm},  with  hm  =  \J\mZm,  m  =  1,2,...  are  orthonormal  basis  in  Hq. 

Then,  using  assumption  (Al)  and  letting  A  =  sup1<n<00  Xn  <  oo,  one  can 

have 

OO  OO 

W(t,u)\1Q  =^2\n\crn(t,Un)\2  <  (  sup  An)  \(Tn{t,Un)\2 

n=l  l<n<oc  n=1 

OO 

<  ^rA^fc2K|2<  A(l  +  ||u||2y), 

n—  1 

which  shows  that  hypothesis  (A. 2)  holds. 

Similarly, 

OO 

\a(t,u)  -  aft,  v)\2Lq  <  A^  \an{t,un)  -  an(t,vn)  |2 

n—  1 

oo 

<  A  K2  ^  kn\Un  ~  Vn\2  =  L\\u  -  v\\y . 

n=  1 

Thus  the  hypothesis  (A. 3.)  holds  true. 

Thus  in  the  abstract  setting  of  stochastic  GOY  model,  the  assumptions  (a.l- 
a.2)  are  required  on  the  noise  coefficient  to  impose  the  hypotheses  (A.1-A.3)  in 
the  Hilbert  space  valued  construction. 

In  the  following  lemma,  we  will  show  that  sum  of  the  linear  and  nonlinear 
operator  is  locally  monotone  in  the  Z4-ball. 

Lemma  3.6.  For  a  given  r  >  0,  let  us  denote  by  B?,  the  closed  l4-ball  in  V : 

®r  =  {v  £  V;  ||u||i4  <  r}  . 

Define  the  nonlinear  operator  F  on  V  by  F(u)  :=  — vAu  —  B{u,u).  Then  for  any 
0  <  £  <  ffp-  where  L  is  the  positive  constant  that  appears  in  the  condition  (A. 3), 
the  pair  (F,  ^ecr)  is  monotone  in  Br,  i.e.  for  any  u  £  V  and  v  £  Br 

r4  „ 

(F(u)  -  F(v),w)  -  ~^\w\2  +  e\a(t,u)  -  a{t,v)\2LQ  <  0, 
where  w  =  u  —  v. 


(3.13) 
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Proof.  First  note  that, 

v(Aw,  w)  =  z'H'U’H2- 

Next  using  the  Lemma  3.4  and  Eq.  (3.10)  from  Lemma  3.3,  we  have 

( B(u,u )  —  B(v,v),w)  =  ( B(v,w )  +  B(w,v)  +  B(w,w),w)  =  (B(w,v),w). 

Now  using  the  definition  of  the  operator  B  and  Eq.  (3.4)  from  Lemma  3.2,  we  get 
for  C  >  0, 


\(B(w,v),w)\  = 


T,  ikr, 

n=  1 

1  *  * 


un-lUJn+lLun 


- 


’n+lVn+2  ' 


n+2Vn+l)wn 


<  c\\v\\i4w\\i4w 

IK 

27 


<  ikikm1/2ikii3/2 


V  ,, 

<2^ 


Since  v  £ 


32  v3 

r,  the  above  relation  yields 


H2IMI  %■ 


-(B{w,v),w)  <  -||w||2 


Hence  by  the  definition  of  the  operator  F, 


(F(u)-F(v),w)<--\\wf  +  ^\w\ 


w„ 


(3.14) 

Finally,  using  condition  (A. 3)  and  that  e  <  ^  we  get  the  desired  result.  □ 

Remark  3.7.  Notice  that  the  trilinear  operator  can  also  be  estimated  as  follows, 


\(B(w,v),w)\  <  CMM 


sup  kn\wn 

l<n<oo 


<  Cf|'u||rc|||'ic||y, 


as  it  is  obvious  that, 

IMIw'1.°°  =  II A1/2tu||00  <  \\A1/2w\\2  =  \\w\\v. 


Then 

\{B(w,v),w)\  <  V- |MI 2  +  ^M2M2. 

Proceeding  as  in  the  previous  Lemma  one  can  see  that  the  operator  F[u)  = 
— uAu  —  B(u ,  u)  is  also  locally  monotone  in  Z2-ball.  But  in  this  paper  we  prefer  to 
work  with  Z4-local  monotonicity,  since  it  is  better  as  the  norm  is  weaker. 
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3.3.  Energy  estimate  and  existence  theory 

Let  Hn  :=  span  {ei,  e2, . . .  ,en}  where  {e^}  is  any  fixed  orthonormal  basis  in  H 
with  each  e,j  £  D(A).  Let  Pn  denote  the  orthogonal  projection  of  H  to  Hn.  Define 
un  =  Pnu ,  not  to  cause  any  confusion  in  notation  with  earlier  un.  Let  Wn  =  PnW. 
Let  <7n  =  Pna.  Define  un'e  as  the  solution  of  the  following  stochastic  differential 
equation  in  the  variational  form  such  that  for  each  v  £  Hn , 

d (un’e(t),v)  =  ( F(un’£(t )),  v)dt  +  (f(t),v)dt  +  y/e(an(t,  un’e(t))dWn(t),v), 

(3.15) 

with  Mn,e(0)  =  Pnu( 0). 

The  standard  theory  of  finite-dimensional  stochastic  differential  equations 
[21]  guarantees  the  existence  of  a  unique  solution  to  (3.15)  under  the  assumptions 
(A.1-A.3)  if  f  £  L2([0,T],H)  and  u(0)  is  Pq  measurable  and  in  L2(P). 


Theorem  3.8.  Under  the  above  mathematical  setting  let  f  be  in  L2([0,  T\,  H),  u(0) 
be  Pq  measurable  and  A|m(0)|2  <  oo.  Let  un,s  denote  the  unique  strong  solution 
of  the  stochastic  differential  equation  (3.15)  in  C([0,T],  Hn).  Then  with  K  as  in 
condition  (A. 2),  the  following  estimates  hold: 

For  all  e  <  27^ ,  and  0  <  t  <  T, 

E\un’*(i)\2  +"  f  E\\u^(s)\\2ds<E\um2  +  -  f  \\f{s)\\2v,ds  +  eKT, 
z  Jo  v  Jo 

(3.16) 

and  for  all  e  <  ^  A  +  \  - 


E 


sup  \un’e{t)\2+  V-  f  \\un*{t)\\2dt  <c(E\um2,  fT  \\m\\v'dt,iy,T)  . 
o <t<T  4  J0  \  Jo  / 


(3.17) 


Also,  for  any  6  >  0  and  e  < 


3  is 
2  K’ 


E\un'e{t)\2e-5t  +  l-  jT  E\\un’E{t)\\2e-5tdt  <  £|u(0)|2 

2  Jo 

+15£m\2e-stdt+£f.  (3.i8) 


Moreover,  if  we  suppose  that  f  £  L4([0, T],H)  and  A|u(0)|4  <  oo,  then  for  all 
e  <  3^  and  0  <  t  <  T, 


E 


sup  \un’e(t)\4e~st  +  Cv  f T  [  |  u" ,£  (t)  1 1 2 1  u" ,£  (t)  | 2  e_  54  dt 
<t<T  Jo 

<  E\u( 0)|4  +  C^T  [T  \\f(t)\\y,e~Stdt  + 

Jo 


eM 

X' 


(3.19) 
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Proof.  Replacing  v  with  un,e  in  (3.15)  and  using  the  properties  of  the  operators 
A  and  B ,  we  notice  that, 

d\un’E(t)\2  +  2is\\un’e(t)\\2dt 

=  2(f(t),un’e(t))dt.  +  £Tr(an(t,un’e(t))Q<rn(t,un’e(t)))dt 
+2V/£(crn(f,  Un’E(t)),Un’e(t))dWn(t). 

Using  the  inequality 

2 ab  <  Sa2  +  ^zb2 
o 

on  2 and  using  the  condition  (A. 2),  we  obtain 

d\u^{t)\2  +  2HK’e(t)||2df  <  HK’£W II2  +  l\\m\\v')dt 
+eK(l  +  \\un’e  (t)\\2)dt  +  2y/e(an(t,un’e(t)),un’s(t))dWn(t). 


Define 


rjv  =  inf|t:  \un’e(t)\2  +  K’e(s)||2ds  >  TvJ  . 

Then  integrating  one  can  deduce 

rtArN 

\un,e(t  A  tjv)|2  +  {v  —  eK)  /  ||Wn>e(s)||2ds 

Jo 

1  rtArN  ftATN 

<H0)\2+-  \\f(s)\\2v,ds  +  eK  ds 

v  Jo  Jo 

rtATN 

+2y/e  /  (crn(s,wn’£(s)),un’£(s))dWn(s). 

Jo 


Let  e  <  Then  taking  expectation  on  both  sides,  and  using  the  fact  that  the 
the  stochastic  integral  appeared  in  the  last  term  of  right  hand  side  in  the  above 
estimate  is  a  martingale,  and  hence  has  a  zero  average,  we  get 


E 


/■tArjv 

^(t/\TN)\2  +  (u-sK)  /  |K'e(s)||2ds 

Jo 


0 

1  rt/\TN 

<E\u(0)\2  +  -  \\f(s)\\2v,ds  +  £KT, 

"  Jo 


thus  we  have  the  result  (3.16). 
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To  prove  (3.17),  we  proceed  in  the  similar  way  as  above,  but  we  take  supre- 
mum  upto  time  T  A  tjv  before  taking  the  expectation, 


E 


/•TAtn 

sup  \un’e(t)\2  +  (v-eK)  /  \\un’e(t)\\2dt 

:<TAtn  Jo 

1 1  pTATN 

<E\u( 0)|2  +  -  /  ||/(i)|&,di 

v  Jo 


0  <t<T  Atn 


+  2  yflE 


sup 

0<t<T  A  rjv 


eKT 


K(S,un-e(S)),Un>E(S))dW„(S) 


(3.20) 


By  means  of  Burkholder-Davis-Gundy  inequality,  condition  (A. 2)  and  Cauchy- 
Schwartz  inequality, 

rt(«7n(s,un’e(s)),un’e(s))dlTn(S) 


sup 

0<t<T  A  rjv 


2  sfeE 

<  2V2 eKE 


(■TAtn  \  1/2 

lo  (i+\\u^m2)\u^(t)\2dt\ 


<  2 V2eKE 


sup  |u 

0<t<T  Arw 


/  pTATN  \ 

;Wl(/  (i  +  IK'E(t)||2)dtj 


1/2' 


< 


(\  pTAtn 

sup  K’E(f)|2  +  £  /  |K’e(f)||2dt  +  T 

0<t<TArN  J  Jo 


■  (3.21) 


Replace  (3.21)  in  (3.20), 


(1  -  s/2eK)E 


sup  \un’e(t)  |2 
0<t<T  Atn 
TAtn 


< 


r 

+  {v-eK-V2eK)  / 
Jo 

1  /*T  ATJV  _ 

£|u(0)|2  +  -  /  ||/(t)||^,dt  + ifT(e+  v^e). 

^  Jo 


Let 


|u"-e(i)||2di 

(3.22) 


e  < 


2R2 


ai./H  +  I-A 


K  2 


Note  TAtn  — >■  T  a.s.  as  AT  — >  00.  Thus  taking  the  limit  in  the  above  estimate 
(3.22)  as  A-  — >  00,  one  can  get  the  desired  energy  estimate  (3.17). 

Next,  we  consider  the  function  e~St[un,e(t)\2  for  S  >  0  and  apply  the  Ito 
Lemma  to  get, 

d  [|u”’E(t) |2e_l5t]  +  2v\\un’e(t) \\2e~5tdt  +  5\un’e(t) \2e~Stdt 
=  [2 (f(t),un’e{t))  +  £Tr(an(t,un’e(t))Qcrn(t,un’E(t)))]  e~Std t 

+2 y/e((Tn(t,  un’e(t)),un’E(t))e~StdWn(t).  (3.23) 
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Note  that 

2(m,u^(t))<5\un’%t)\2+lsm\2. 

Hence  upon  writing  (3.23)  in  the  integral  form,  taking  expectation  and  using 
condition  (A. 2),  one  can  get 

E\un’E{i)\2e-5t  +  2vE  [  \\un,e  (t)\\2  e~St  dt 

Jo 

<E\u(0)\2+l  f  \f(t)\2e-5tdt  +  eKE  [  (1  +  \\un’E(t)\\2e~Stdt, 

°  Jo  Jo 

which  yields  the  estimate  (3.18)  for  all  e  < 

To  prove  (3.19),  we  first  apply  Ito  Lemma  on  the  function  |itn’e(t)|4e-,5t  to 
get, 

d[|un’e(f)|4e-5t]  +  4i'||Mra’e(t)||2|Mn’£(<)|2e-<5tdf  +  <5|'it7l,e(f)|4e-'5tdt 

=  \un’£(t)\2[4(f(t),un^(t))+8eTv(an(t,un^(t))Qan(t,un’E(t)))]e-5tdt 
+^(an(t,un'£(t)),un’e(t))\un’e(t)\2e-5tdWn(t).  (3.24) 

Using  the  fact  that 

4(/(t),U”’£(t))K’£(t)]2  <  Cs\f(t)\4  +  6\un’E(t)\\ 
and  applying  the  condition  (A. 2)  and  integrating  we  have, 

\un'E(t)\Ae~5t +  (Av-SeK)  [  ||iin'e(a)||2|un>e(«)|2e-4ads 

Jo 

<  E\u(0)\4  +  Cs  f  \f(s)\Ae-Ssds 
Jo 

+4 Ve  [  (an(s,un’e(s)),un’s(s))\un’s(s)\2e-SsdWn(s). 

Jo 

Finally,  taking  supremum  in  0  <  t  <  T,  then  taking  expectation  on  both  sides  and 
using  the  Burkholder-Davis-Gundy  inequality  on  the  stochastic  integral  term,  we 
get  the  estimate  (3.19).  □ 

Definition  3.9.  (Strong  Solution)  A  strong  solution  u£  is  defined  on  a  given  proba¬ 
bility  space  (fl,E,Et,P)  as  a  L2(U;  L°°(0, T;  H)nL2(0,T;  F)nC(0,T;  H))  valued 
adapted  process  which  satisfies  the  stochastic  GOY  model 

du£  +  [vAue  +  B(ue,  ue )]  dt  =  f(t)dt  +  y/ eo(t ,  u,e)dW{t)  (3.25) 

me(0)  =  u0, 

in  the  weak  sense  and  also  the  energy  inequalities  in  Theorem  3.8. 

Monotonicty  arguments  were  first  used  by  Krylov  and  Rozovskii  [22]  to  prove 
the  existence  and  uniqueness  of  the  strong  solutions  for  a  wide  class  of  stochastic 
evolution  equations  (under  certain  assumptions  on  the  drift  and  diffusion  coeffi¬ 
cients),  which  in  fact  is  the  refinement  of  the  previous  results  by  Pardoux  [30,31] 
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and  also  the  generalization  of  the  results  by  Bensoussan  and  Temam  [3] .  Menaldi 
and  Sritharan  [27]  further  developed  this  theory  for  the  case  when  the  sum  of  the 
linear  and  nonlinear  operators  are  locally  monotone. 

Theorem  3.10.  Let  u(0)  be  Tq  measurable  and  F|uo|4  <  oo.  Let  f  £  L4(0,T;  V'). 
We  also  assume  that  0  <  e  <  f  and  the  diffusion  coefficient  satisfies  the  condi¬ 
tions  (A.l  A. 3).  Then  there  exists  a  unique  adapted  process  uE(t,x,w)  with  the 
regularity 


ue  £  L2(0;C(0,T;l?)nL2(0,T;  V)) 

satisfying  the  stochastic  GOY  model  (3.25)  and  the  a  priori  bounds  in 
Theorem  3.8. 

Proof.  Part  I.  (Existence) 

Using  the  a  priori  estimate  in  the  Theorem  3.8,  it  follows  from  the  Banach-Alao- 
glu  theorem  that  along  a  subsequence,  the  Galerkin  approximations  {un,E}  have 
the  following  limits: 

un,e  — >  ue  weak  star  in  L4(U;  L°°(0 ,T;H))  fl  L2(fl;  L2(0,  T;  V)), 
F{un’E) — *  Fq  weakly  inL2(f];  L2(0,  T;  V')),  (3.26) 

<Jn{-,un'e)  —  weakly  in  L2(0;  L2(0,  T;  Lq)). 

The  assertion  of  the  second  statement  holds  since  F(un,E)  is  bounded  in  L2(fl; 
L2(0,  T;  U')).  Likewise  since  diffusion  coefficient  has  the  linear  growth  property 
and  un,e  is  bounded  in  L2(0,  T;  V)  uniformly  in  n,  the  last  statement  holds.  Then 
uE  has  the  Ito  differential 

dne(t)  =  Fq  (f)d£  +  y/eSe(t)dW(t)  weakly  in  L2(U;  L2(0,  T;  V')). 

Let  us  set, 


r(f):=4  f  \\vs{s)\\Us,  (3.27) 

v  Jo 

where  ve(t,x,uj)  is  any  adapted  process  in  L°°(fl  x  (0 ,T);H).  Here  we  suppress 
the  dependence  of  e  in  the  notation  of  r  to  make  it  easier  to  read.  Then  applying 
the  Ito  Lemma  to  the  function  2e_r^^|Mrl’e(t)|2,  one  obtains 


d 


i tr(t)K’e(f)]2 


=  e-r(t)  (2F(un-e(t))  -r(f)un’e{t),un’E(t))dt 
+£e~r{t)\an(t,un’E(t))\lQdt 
+2^e~rW  ( an{t ,  un’e(t)),  un’E(t))  dW(t). 
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Integrating  between  0  <  t  <  T  and  taking  expectation, 


E 


e~r^\un’E(T)\2  -  \un’e(0)\2 

T 


=  E 


[  e~r{t)  (2 F(un’E(t))  -  f(t)un’e(t),un'e(t))  d t 

Jo 

+eE  [  e~rW\an{t:un’s(t))\l  dt 
Jo 

+2y/eE  [  e-rW(an(t,un'E(t)),un’e(t))dW(t). 

Jo 


The  last  term  on  the  right  hand  side  vanishes  since  the  integral  inside  the 
expectation  is  a  martingale.  Then  by  the  lower  semi-continuity  property  of  the 
weak  convergence, 


lim  inf  E 


e~rW  (2F(un’E(t))  -f(t)un’E(t),un’e(t))dt 


+e  [  e  rW\crn(t,un'E{t))\l  dt 
Jo 


=  lim  inf  E 


e~riT)\un'E{T)\2  -  \un’E{0)\2 


>  E 

=  E 


e~r^\uE(T)\2  -\uE(0)\2 


=-r(t) 


(2-FoW  -r(t)uE(t),uE{t))dt  +  e  [  e  rW\SE\l  dt 

Jo 


(3.28) 


Now  by  monotonicity  property  from  Lemma  3.6, 


2  E 


f  e"r(t)  ( F(un'E(t ))  -  F(vE(t)),  un'e(t)  -  vE{t))dt 
Jo 


-E 

+eE 
<  0. 


e-rWf(t)\un’E(t)-vE{t)\2dt 


f  e  r(t)\an{t,un’E{t))  -  an{t,ve(t))\lQdt 
Jo 
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Rearranging  the  terms, 

f 

la 


E 


<  E 


e_r(*)  (2 F(un’E(t))  -  f(t)un’e(t),un'E(t))dt 

+£  [  err(t)\an(t,un’e{t))\lQdt 
Jo 

rT 

/  e~r(t)  (2 F(un’e{t))  -  r(t)(2un'e(t)  -  ve(t)),  ve(t))  d t 

Jo 


+E 


+  eE 


e-rW  (2F(ve(t)),un’e{t)-v£(t))dt 


J  e  r(t)  (2an{t,un'£{t))  -  an(t,vs(t)),an(t,ve(t)))LQ  dt 

Taking  limit  in  n,  using  the  result  from  (3.28)  and  rearranging,  we  obtain 

’  rT 

E  /  e~r^  (2F^(t)-2F(v£(t)),u%t)-ve(t))dt 
Jo 


+E 

+eE 
<  0. 


e~rWf(t)\ue(t)  ~v£(t)\2dt 


e-r^\\S(t)-a(t,v^t))\\t0dt 


Notice  that  for  ve  =  ue ,  S(t)  =  a(t,ue(t)).  Take  ve  =  ue  —  A we  with  A  >  0  and 
we  is  an  adapted  process  in  L2(fi;  C(0,T;  H)  flL2(0,T;  V))  Then, 


A  E 


o~r{t)  {2F^t)-2F(ue-\we)(t),we{t))dt  +  \  f  e~r{t)f(t)\we{t)\2dt 

Jo 


<  0. 


Dividing  by  A  on  both  sides  of  the  inequality  above  and  letting  A  go  to  0,  one 
obtains 


E 


(Fo(t)  -  F(uE{t)),we(t))dt 


<  0. 


Since  we  is  arbitrary,  we  conclude  that  Fg(f)  =  F(ue(t)).  Thus  the  existence  of 
the  strong  solution  of  the  stochastic  GOY  model  (3.25)  has  been  proved. 
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Part  II.  (Uniqueness) 

If  ve  £  L2(fl;  C(0,T;H)  n  L2(0,  T;  V))  be  another  solution  of  the  equation  (3.25) 
then  wE  =  ue  —  ve  solves  the  stochastic  differential  equation  in  L2(fl;  L2(0,  T;  V')), 

du»E(f)  =  (F(uE(t))  —  F(vE(t)))dt  +  \Je(a(t,uE{t))  —  <r(t,  vE(t)))dW(t).  (3.29) 

We  denote  ad  =  a(t,uE(t))  —  a(t,vE(t)). 

We  now  apply  Ito  Lemma  to  the  function  2e~r^t') \wE (t)\2  and  using  the  local 
monotonicity  of  the  sum  of  the  linear  and  nonlinear  operators  A  and  B ,  e.g. 
equation  (3.14),  we  get 

or4 

e-rWd\wE(t)\2+ve-rW\\wE(t)\\2dt  <  -^e-r{t)\wE(t)\2dt+ee-rWTiiadQad)dt 

+2y/ee-rW(ad,wE{t))dW{t).  (3.30) 

Using  condition  (A. 3), 

d(e-r{t)|toe(i)|2)  +  ve~rW\\wE(t)\\2dt  <  eLe-r^\\wE{t)\\2dt 

+2y/e :e~r('t\ad,  wE (t))dW (t) . 

(3.31) 


Finally  integrating  in  0  <  t  <  T,  taking  expectation  on  both  sides  and  noting 
£  <  jr  and  the  fact  that 


D-r(t) 


(ad,we(t))dW(t) 


is  a  martingale  for  T  <  oo,  we  obtain  P-a.s. 


E 


>-rW\wE(t)\2  <E\w(0)\2, 


which  assures  the  uniqueness  of  the  strong  solution. 


□ 


4.  Large  deviation  principle  continued 

Let  us  recall  the  stochastic  GOY  model  in  consideration, 

ditE  +  [vAue  +  B(uE ,  uE)}  dt  =  f(t)dt  +  y/ea(t,  uE)dW(t)  (4.1) 

«E(0)  =  £, 

has  a  unique  strong  solution  in  the  Polish  space  X  =  C([0,T];  iJ)  n  L2(0,T;  V). 
The  solution  to  the  stochastic  GOY  model,  denoted  by  uE  can  be  written  as 
Ge{W(-)),  for  a  Borel  measurable  function  QE  :  C'([0,T];ff)  — >  X  (see  [21],  page 
310,  [38],  Chapter  X,  Corollary  4.2).  The  aim  of  this  section  is  to  verify  that  such 
a  QE  satisfies  Assumption  2.10.  Then  applying  the  Theorem  2.11  the  LDP  for 
{uE  :  e  >  0}  in  X  can  be  established. 

The  LDP  for  {ue  :  e  >  0}in  X  have  been  proved  here  systematically  in  four 
steps.  In  the  first  and  second  Theorems  we  show  the  well  posedness  of  certain 
controlled  stochastic  and  controlled  deterministic  equations  in  X.  These  results 
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help  to  prove  the  last  two  main  Theorems  on  the  compactness  of  the  level  sets  and 
weak  convergence  of  the  stochastic  control  equation  stated  in  Assumption  2.10. 


Theorem  4.1.  For  any  v  £  Am,  0  <  M  <  oo,  the  stochastic  control  equation 
d  u%(t)  +  [vAu%{t)  +  B{uev(t),u%{t))\dt  =  [f(t)  +  d{t,uev(t))v{t)]dt 

+Vea(t,uEv(t))dW(t),  (4.2) 

<(0)  =  eeir, 

has  a  unique  strong  solution  in  L2(fl;X),  where  X  =  C(0,T;H)  flL2(0,T;  k), 
/  £  L4(0,T;  V')  and  a,  a  both  satisfy  the  hypotheses  A.l-  A. 3.  in  Sect.  3. 

Proof.  We  first  prove  that  if  uEv{t)  is  a  strong  solution  of  the  stochastic  controlled 
equation  (4.2),  the  following  energy  estimate  holds: 

e(  sup  K(t)|2+  [  \\uev(t)\\2dt  )  <  C,  (4.3) 

\0<t<T  Jo  J 

where 

is  a  positive  constant. 

To  prove  the  above  estimate,  we  apply  Ito  formula  on  \uev(t)\2  and  integrating 
in  0  <  t  <  T, 

K(t)\2  +  u  f  \\ul(S)\\2ds 
Jo 

<IC|2  +  -/  \\f(s)\\v'ds  +  2  [  (d(sX(s)Ms)X(5))  ds 

v  Jo  Jo 

+  £  f  Tr(a(s,<(s))Qcr(s,<(s)))ds+2Ve  [  (a(s,uev(s)),uev(s))  dW(s). 
Jo  Jo 

(4.4) 

Notice  that, 

2  [  (d(sX(s)Ms)X(s))ds 
Jo 

<2  f  |CT(sX(s))|LQKs)|oK(s)|ds 
Jo 

~i^JQ  I^OXO))!2 ods+x/  Ks)loK(s)l2ds 

, ,  P  AK  P 

<-  /  (1+X(s)||2)ds+—  /  |X)loKWI2ds 
^  Jo  v  Jo 

<j[  (1+ IKMlDd s+AI^AJ  sup  K(i) |2,  (4.5) 

4  Jo  V  0 <t<T 
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and  for  e  <  -At, 


4  K’ 

rt  rt 

£  Jo  Tr(a(s’uv(s^®a(s’uv(s^ds  -  4  /  (1  +  IK0)ll2)ds-  (4-6) 

After  rearrangement  of  the  Eq.  (4.4),  we  take  supremum  in  time  0  <  t  <  T  and 
then  expectation  to  get, 

rT 


E 


sup  |<(t)  | 
0<t<T 


V 

2 


m\\2dt 


<icis 


z/  1 
2T+u 


wmwv’dt- 


4KM 


-E 


sup  \uev{t)Y 

.0  <t<T 


+  2  yfiE 


sup 

o  <t<T 


(cr(s,<(s)),<(s))dTE(s) 


(4.7) 


The  last  term  of  the  above  equation  can  be  estimated  in  a  similar  manner  as  in 
(3.21), 

rt 


2  sfeE 


sup 

0<£<T 


(cr(s,<(s)),<(s))dlE(s) 


<  V2eK 


E(  sup  \uev{t)\2)  +  E  [T  \\uim2 At  +  T 
0  <t<T  Jo 


(4.8) 


Replacing  (4.8)  in  (4.7),  and  considering 

.  f  v  v 2  (1  -  AKM/v)'2 

nil  <  -  -  - 


e  <  mm  • 


2  K2 


(4 K1  8K2  ’ 
we  get  the  energy  estimate  (4.3). 

The  proof  of  the  existence  and  uniqueness  of  the  strong  solution  of  the 
stochastic  controlled  equation  (4.2)  follow  from  the  Theorem  3.10,  only  a  few 
modifications  are  needed  due  to  the  presence  of  the  control  term.  □ 

Theorem  4.2.  Let  v  £  L2(0,  T;  Ho),  f  £  L4(0  ,T\V')  and,  a  satisfy  the  hypotheses 
A.1.-A.3.  in  Sect.  3.  Then  the  equation 

d uv(t)  +  \vAuv(t)  +  B(uv(t),uv(t))]dt  =  f(t)dt  +  a(t,  uv(t))v(t)dt,  (4.9) 
where  uv(fS)  =  £  £  H,  has  a  unique  strong  solution  in  X  =  C(0,  T;  i7)flL2(0,  T;  V). 

Proof.  This  result  can  be  considered  as  a  particular  case  of  the  previous  Theorem 
4.1,  where  the  diffusion  coefficient  is  absent.  □ 


Next  we  state  a  important  lemma  from  Budhiraja  and  Dupuis  [4]. 

Lemma  4.3.  Let  {vn}  be  a  sequence  of  elements  from  Am  for  some  finite  M  > 
0.  Let  vn  converges  in  distribution  to  v  with  respect  to  the  weak  topology  on 
L2(0,T;77o).  Then  fQ  vn(s)ds  converges  in  distribution  as  C(0,T;  H)-valued  pro¬ 
cesses  to  fQ  v(s)ds  as  n  —>  oo. 

Now  we  are  ready  to  check  the  Assumptions  2.10. 
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Theorem  4.4  (Compactness).  Let  M  <  oo  be  a  fixed  positive  number.  Let 

Km  :={uv£  C(0,  T;  H)  fl  L2 (0,  T;  V);  v  £  SM  } , 

where  uv  is  the  unique  solution  in  X  =  C(0,T-,H)  nL2(0,T;’C)  of  the  deter¬ 
ministic  controlled  equation  (4.9),  with  uv{ 0)  =  £  £  H .  Then  Km  is  compact 
in  X . 


Proof.  Let  us  consider  a  sequence  {uVn}  in  Km,  where  uVfl  corresponds  to  the 
solution  of  (4.9)  with  control  vn  £  Sm  in  place  of  v,  i.e. 

d uVn  (t)  +  \vAuVn(t)  +  B(uVn  (t),  u„n(f))]df  =  f(t)dt  +  oft,  uVn  (t))vn(t)dt,  (4.10) 

with  u„n(0)  =  £  £  H.  Then  by  weak  compactness  of  Sm,  there  exists  a  sub¬ 
sequence  of  {vn},  still  denoted  by  {un},  which  converges  weakly  to  v  £  Sm  in 
L2(0  ,T;H0). 

We  need  to  prove  uVn  — >  uv  in  X  as  n  — >  oo,  or  in  other  words, 

sup  \uVn(t)  -  uv(t)\2  +  I  ||u„n(t)  -  uv(t)\\2dt  ->  0,  (4.11) 

o  <t<T  Jo 

as  n  — >  oo. 

According  to  the  Theorem  4.2,  uv  is  unique  strong  solution  in  X  of  the  deter¬ 
ministic  controlled  equation  (4.9).  Hence  it  is  obvious  to  note  that,  uv  satisfies 
the  following  a-priori  estimate 

sup  \uv(t)\2  +  [  \\uv{t)\\2dt  <  C,  (4.12) 

0  <t<T  Jo 

where 

C  =  C  (|£|2,  \\f\\v'dt,  v,  K,T,M^J 

is  a  positive  constant. 

For  the  proof,  we  refer  the  Theorem  4.1,  where  the  stochastic  version  of  the 
above  a  priori  estimate  has  been  worked  out. 

Let  wVn  =  uVn  —  uv.  Then  wVn  satisfies  the  following  differential  equation 

dwVn(t)  +  [vAwVn(t)  +  B(uVn(t),uVn(t))  -  B(uv(t),uv(t))]dt 

=  [o(t,uVn(t))vn(t)  -  o(t,uv(t))v(t)]dt,  (4.13) 


which  yields 


vn(t) |2  +  2^  [  ||w„„(s)||2ds 

Jo 

+2  /  (B(uVn(s),uVn(s))  -  B(uv(s),uv(s)),wVn(s))ds 
Jo 

=  2  /  (o(s,uVn(s))vn(s)  -  o(s,uv(s))v(s),wVn(s))ds. 

Jo 


(4.14) 
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First  note  that,  from  Lemma  3.4, 

B(^^vn  ’  )  B(\LV ,  'Uv')  —  33 (uv ,  WVrl )  “b  B(wVn  ,  Uv  )  “b  B(wVn  ,  'Wvn  )  ■ 

Using  the  above  expression  and  the  properties  ( ii )  and  ( iv )  of  the  bilinear  operator 
B  given  in  Lemma  3.3,  one  can  find 

2  \(B(uVn(s),uVn(s))  -  B(uv(s),uv(s)),wVn(s))\ 

=  2  \  (B(wVn(s),  uv(s)),wVn  (s))  | 

<  2||S(u;Wn(s),u„(s))||1v'|kvn(s)|| 

<2|tu„n(s)|  |u„(s)|  ||w„„(s)|| 

<  \\ K„(s)||2  +  ^K„(s)|2K(s)|2-  (4.15) 


Also  notice  that, 


(a(s,uVn(s))vn(s)  -  <r(s,  uv(s))v(s) ,  wVn(s))  ds 


<  /  |((cr(s,u„n(s))  -  cr(s,ut,(s)))  vn(s),  W„n(s))|  ds 
Jo 


(a(s,uv(s))(vn(s)  -  w(s)),u;t,rl(s))ds 


<Vl  ||ttb;n(s)||  |u>„n(s)]  |v„(s)|0ds 
Jo 


(cr(s,  uv(s))(vn(s)  -  w(s)),wh,„(s))ds 


v 

<  - 
“  4 


[  \\wvn{s)\\‘2ds+--  [  |u>„n(s)|2  K(s)|ods 

Jo  v  Jo 


sup 

0  <t<T 


(<r(s,uv(s))(vn(s)  -  v(s)),wVn(s))ds 


(4.16) 


By  the  boundedness  of  {|wh,„(s)|2}  in  C(0,T;17),  and  using  the  Lemma  4.3,  the 
second  integral  on  the  right  side  of  (4.16)  goes  to  0  as  n  — >  oo.  Therefore,  given 
any  e  >  0,  there  exists  an  integer  N  large  so  that  for  all  n>  N, 


sup 

0<t<T 


(a(s,uv(s))(vn(s)  -  v(s)),ot)„„(s))  ds 


<  e/2. 


(4.17) 


Consider, 


C]jV  =  max 


2  2  L 


v  v 
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Applying  (4.15),  (4.16)  and  (4.17)  in  (4.14),  one  obtains  for  n>  N, 

\wvn(t)\2  +  v  f  ||«An(s)||2ds  <  CL,V  [  l^„(s)|2(|u^(s)|2  +  |un(s)|^)  da- 
Jo  Jo 


(4.18) 


Hence  by  Gronwall’s  inequality, 

rT 


sup  \wVn(t)\2  +  v  f  ||w^W||2dt<  eeC£-^oT(|^W|2+l^(t)|g+1)dt_  (410) 
0  <t<T  J  0 

The  arbitrariness  of  e  finishes  the  proof.  □ 

Remark  4.5.  From  Theorem  4.1  one  can  see  that  the  equation 

du%s{t)  +  [Au%e{t)  +  B(ule(t),ule(t))]dt 

=  [f{t)  +  a(t,ueve(t))v£(t))dt  +  y/ea(t,ueve{t))dW{t),  (4.20) 

with  rtEe(0)  =  £  £  H,  has  unique  strong  solution  in  L2(fl;  X). 

As  before,  the  solution  admits  a  representation  usve  =  Ge(W(-)  + 

fQ  ve(s)ds)  by  pathwise  uniqueness  of  the  solution,  and  the  Girsanov  theorem. 

For  all  v  £  L2(0,T;  Hq),  let  uv  be  the  solution  of  the  deterministic  control 
equation 

d  uv(t)  +  \vAuv(t)  +  B(uv(t),uv(t))]dt  =  f(t)dt  +  o(t ,  uv(t))v(t)dt, 

with  initial  condition  u„(0)  —  ^  G  H. 

Note  that  fQ  v(s)ds  £  C([0,T]  :  H0).  Define  G°  :  C([0,T]  :  H0)  ->  X  by 

G°(h)  =  uv  if  h  =  j  w(s)ds  for  some  v  £  L2([0,T]  :  H0). 

Jo 

If  h  cannot  be  represented  as  above,  then  define  G°(h)  =  0. 

Theorem  4.6  (Weak  convergence).  Let  {wE  :  £>0}  C  Am  converges  in 
distribution  to  v  with  respect  to  the  weak  topology  on  L2(0,  T;  Hq)  .  Then  GE(W(-)+ 
-A  fg  ve(s)ds )  converges  in  distribution  to  G°(f0  v(s)ds)  in  X,  as  e  — >•  0. 

Proof.  Since  Sm  is  Polish,  the  Skorokhod  representation  theorem  can  be  intro¬ 
duced  to  construct  processes  (v£,v,We)  such  that  the  distribution  of  (v£,v,  W£) 
is  same  as  that  of  (v£,v,W),  and  v£  — >  v  a.s.  in  the  weak  topology  of  Sm-  Thus 
/0t£iE(s)ds  ->  foV(s)ds  weakly  in  Hq  a.s.  for  all  t  £  [0,T].  Without  any  loss  of 
generality,  we  will  write  (vE,v,W)  in  what  follows,  though  strictly  speaking,  one 
should  write  (v£,  v,  W£). 

Let  w£s{t )  =  ueve(t)  —  uv(t).  We  need  to  prove, 

sup  \w£s{t)\2  +  f  ||wEe(i)||2dt  ->  0  (4.21) 

0  <t<T  Jo 

in  probability  as  e  — >  0. 
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Notice  that,  applying  similar  estimate  as  in  previous  Theorem,  equation 
(4.18),  one  can  get 


< 


(t)\2  +  v  f  \\wUs)fds 
Jo 

3 CL>V  f  |<e(s)|2  (K(s)|2+  |we(s)|q  +  l)ds 
Jo 

+  [  \a(s,uv(s))(vE(s)  -  v(s))\2ds  +  eK  f  (1  +  ||<e(s)||2)ds 
Jo  Jo 

+2y/e  f  (a{s,ueV'(s)),wl*{s))dW{s). 

Jo 


(4.22) 


We  take  supremum  in  0  <  t  <  T,  then  expectation  on  the  above  inequality,  and 
use  similar  estimate  (with  the  help  of  Burkholder-Davis-Gundy  inequality)  as  in 
(3.21)  on  the  last  term  of  the  right  hand  side  to  get, 


E 


I  sup  \w£ve(t)\2  (K(t)|2  +  \ve{t)\l  +  1)  df 
0  0<t<T 


sup  \wEve(t)  \‘ 
0  <t<T 


sup  \wSve[t)  \2+V  [  \\weve(t)\\2dt 
0<t<T  Jo 

<  ZCl,vE 

+  I  \cr(t,uv(t))(vE(t)  -  v(t))\2dt  + eK(C +  T) 

Jo 

+V2eK  (c  +  T+E 

Assume  that  e  <  97^2  •  Then  the  Gronwall  inequality  yields 

sup  \wSve{t)\2  +  V  [  \\wEve{t)\\2dt 

)<t<T  Jo 

<  {(e+V2e)K{C  +  T)  +  J  \a(t,uv(t))(ve(t)  —  v(t))\2dt 


E 


xe 


,3 CL,„Jo  (|ii„(i)|2  +  hep)lo  +  1)dt 


(4.23) 


(4.24) 


Since  vE  — >  v  a.s.  in  the  weak  topology  of  Sm,  it  is  clear  from  the  Eq.  (4.24)  that 
as  e  — >  0, 


E 


sup  \wEve(t)\2 +  v  (  \\wEvs(t)\\2dt 
0  <t<T  Jo 


0. 


518 


U.  Manna  et  al. 


NoDEA 


Let  S  >  0  be  any  arbitrary  number.  Then  by  Markov’s  inequality 


0  as  e  — >  0. 


o  <t<T  Jo 

in  probability  as  e  — >  0.  The  proof  is  now  complete. 


□ 


Remark  4.7.  Sabra  shell  model  of  turbulence  is  the  other  well  accepted  model  in 
the  literature,  and  the  fundamental  difference  with  the  GOY  model  lies  in  the 
number  of  complex  conjugation  operators  used  in  the  nonlinear  terms  which  are 
responsible  for  differences  in  the  phase  symmetries  of  the  two  models,  and  as  a 
consequence,  Sabra  shell  model  exhibits  shorter-ranged  correlations  than  the  GOY 
model  (see  [24]).  The  equations  of  motion  of  the  stochastic  Sabra  shell  model  have 
the  following  form 


U.  U/yj  .  n  •  /  7  +  ■  ■ 

- -  4-  4-  7,  ( 77,„_i_t77,.-  ,  ,  4-  bk,„u„A_iu, 


'n—  1 


for  n  =  1, 2, . . . , 


along  with  the  boundary  conditions 

U- 1  =  uq  =  0. 


Under  the  same  assumptions  on  the  noise  and  noise  coefficient  given  in  Sect.  3, 
and  under  the  same  functional  setting,  the  existence  and  uniqueness  of  the  strong 
solution  can  be  established  in  L2(fl;  C(0,T;  H)  D  L2(0,  T;  V)).  Moreover,  by  pro¬ 
ceeding  in  the  similar  fashion  as  in  Sect.  4,  one  can  easily  verify  the  key  estimates 
and  prove  the  large  deviation  principle  for  the  solution  of  the  stochastic  Sabra 
model  in  the  Polish  space  (7(0,  T;  H)  flL2(0,T;  V). 
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